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Abstract 

We develop a synthesis of Turing's paradigm of computation and 
von Neumann's quantum logic to serve as a model for quantum com- 
putation with recursion, such that potentially non-terminating com- 
putation can take place, as in a quantum Turing machine. This model 
is based on the extension of von Neumann's quantum logic to par- 
tial states, defined here as sub-probability measures on the Hilbert 
space, equipped with the natural pointwise partial ordering. The sub- 
probability measures allow a certain probability for the non-termination 
of the computation. We then derive an extension of Gleason's theorem 
and show that, for Hilbert spaces of dimension greater than two, the 
partial order of sub-probability measures is order isomorphic with the 
collection of partial density operators, i.e. trace class positive operators 
with trace between zero and one, equipped with the usual partial order- 
ing induced from positive operators. We show that the expected value 
of a bounded observable with respect to a partial state can be defined 
as a closed bounded interval, which extends the classical definition of 
expected value. 

1 Introduction 

The rigorous mathematical foundation of quantum mechanics is generally 
agreed to be based on von Neumann's formulation, which uses the notion of 
a state of quantum logic, a probability measure on the collection of quantum 
events, i.e., on the closed subspaces of the Hilbert space ^2 Ell- An observ- 
able, discrete or continuous, is defined as a mapping from the Borel subsets 
of the real line to the collection of states of the quantum logic. These two 
notions are then used to derive the expected value of an observable, which 
provides a consistent quantum theory treating discrete and continuous ob- 
servables in a uniform manner. Gleason's fundamental theorem shows that, 
for Hilbert spaces of dimension greater than two, the states of logic are in 
one to one correspondence with density operators, i.e. trace class operators 
on the Hilbert space with trace one, and thus allows us to work with the 
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more convenient density operators instead of the probability measures on 
the collection of closed subspaces of the Hilbert space. 

In this paper, we aim to synthesize von Neumann's formulation of quan- 
tum mechanics with Turing's paradigm of computation, in which classical 
iteration or recursion plays a foundational role. 

All standard programming languages allow for non-trivial recursion, a 
fundamental feature of computation which can potentially result in the non- 
termination of programs. In a pioneering paper in 1936, Alan Turing proved 
that it is not decidable in general if a program terminates, i.e. one cannot in 
general determine, in a finite amount of time, if a given program terminates 
on a given input [El El- This extremely important result in computer science 
has since been called "the halting problem". It is a basic consequence of 
the halting problem that in probabilistic programs, when one computes the 
probability of the outcome of an event as in quantum computation, sub- 
probability measures, rather than probability measures, are used to model 
the probabilistic results: the total probabilities of the definite outcome add 
up to a number less than one, allowing a certain probability for the non- 
termination of the program; see for example 0] . 

In order to reconcile von Neumann's quantum logic with the halting 
problem, we therefore consider sub-probability measures on quantum events, 
i.e. measures fj, on the closed subspaces of the Hilbert space TL such that 
< n(TC) < 1, which we call partial states of the quantum logic. There 
is a natural partial order on partial states: /xi < fi2 if fJ'i(A) < A*2(^4) 
for all closed subspace A of 7i, i.e., for any quantum event the probability 
of outcome with respect to Hi is at least as much as that of which 
implies that the probability of non-termination for fi2 is at most that of H\. 
This partial order is complete in that any increasing chain has a supremum, 
more generally it is directed complete in the sense that any directed set has 
a supremum. Recall that a non-empty subset of a partially ordered set is 
directed if for any two elements in the subset there is an element of the subset 
above both; thus a directed subset is a generalization of an increasing chain 
in a partial order. A directed set of partial states represents a consistent 
set of computations: for any pair of partial states in the set, there exists a 
partial state in the set that assigns to each quantum event a probability at 
least as great as either of those assigned to the quantum event by the two 
partial states. 

We then consider trace class positive operators with trace between zero 
and one, which we call partial density operators, equipped with the partial 
ordering given by B < A if A — B is a positive operator. This partial order 
is also directed complete, i.e., any directed subset has a supremum. We 
show that, for Hilbert spaces of dimension greater than two, Gleason's onto 
map from the collection of density operators to that of states of quantum 
logic extends to an order isomorphism between the directed complete partial 
order of partial density operators and that of the partial states of the logic, 
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preserving the supremum of directed subsets. Any directed complete partial 
order is equipped with a natural To topology, the so-called Scott topology, 
which corresponds to the convergence of a directed set to its supremum. The 
extended Gleason map is then continuous with respect to the Scott topolo- 
gies of the directed complete partial orders of the partial density operators 
and the partial states of the quantum logic. All in all, this shows that par- 
tial density operators provide a convenient model for quantum computation 
with recursion. The partial order of partial density operators has already 
been used in ^ to develop a functional programming language for quan- 
tum computation. In a typical recursive computation, the partial density 
operator in each loop of iteration increases in the partial order and in the 
limit one obtains the supremum of this increasing chain of partial density 
operators, which can indeed be a usual density operator; see [51 Section 5.5] 
for an example. 

Finally, we consider any bounded observable, with a discrete or contin- 
uous spectrum, given as a map from the Borel subsets of the real line to 
the lattice of closed subsets of the Hilbert space. We define the expected 
value of the observable with respect to a partial state as a compact real 
interval and show that, for a fixed bounded observable, the expected value 
as a function from the directed complete partial order of partial states to 
the directed complete partial order of the collection of non-empty compact 
intervals ordered by reverse inclusion is Scott continuous. By the spectral 
theory of operators, the bounded observables are in one-to-one correspon- 
dence with bounded self-adjoint operators on the Hilbert space. We thus 
obtain the expected value of a bounded self adjoint operator with respect 
to a partial density matrix and show that for a fixed partial density op- 
erator, the expected value, as an interval-valued map, depends linearly on 
commuting bounded self adjoint operators. 

2 Partial states of quantum logic 

To fix our notations, we briefly review von Neumann's formulation of the 
states of quantum logic. Let TL be a separable (finite or infinite dimensional) 
Hilbert space over the complex numbers. We denote the set of closed sub- 
spaces of Tt partially ordered by inclusion by S(Tt). A closed subspace is 
also referred to as a quantum event. Then S(TC) is a non-distributive com- 
plete lattice with A A B = A n B and A V B = span{j4, B}, the subspace 
generated by A and B. 

A probability measure on S(TC) is a mapping p : S(7i) — > [0, 1] such that 

(i) p({0}) = o, 

(ii) p(H) = 1, 
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(iii) p(Vn>o A n) = J2 n >oP( A n), for any sequence A n of mutually orthogo- 
nal subspaces A n . 

A state of the quantum logic or of quantum events is a probability mea- 
sure on S(H); it gives the probability of each event A £ S(H) and thus all the 
information required to compute any physical property of the system. We 
now consider any quantum programming language, which computes proba- 
bilities of quantum events and is allowed to use recursion. A partial state of 
the quantum logic is defined as a subprobability measure of S(H), i.e., as a 
mapping p : S(H) — ► [0, 1] such that 

(i) p({0}) = 0, 

(ii) p(H) < 1, 

(iii) p(Vn>o A n) = !>2n>oP( A n), for any sequence A n of mutually orthogo- 
nal subspaces. 

The number 1 — p{H) > is the probability of non-termination of the 
program or the computation process. We note that this number, and hence 
the probability of termination, can be non-computable when the Hilbert 
space Ti is infinite dimensional. 

In fact, for an orthonormal basis (ej)j>o of H and any non-computable 
number a between and 1, there is a subprobability measure p on 7i with 
p(Tt) = a, such that p assigns a rational number to any finite subspace 
generated by (ej)j>o. Such p can be defined as follows. We consider the 
binary expansion of a: 

i>0 

where € {0, 1} for i > O.and define the partial state p by its values on the 
basis vectors as p(span{ej}) = ^fr- Then p assigns a rational number to 
any finite dimensional subspace generated by the orthonormal basis (ej)j>o 
and we have p(H) = a. 

There is a natural notion of partial order on partial states. For two 
partial states, we define p\ < P2, if for all quantum events A G S(7i), we 
have p\{A) < P2^A), i.e., p2 gives more probability to any quantum event 
and is more likely to terminate compared to p\. We need the following 
lemma, which is a generalization of a corresponding well-known result for 
double sequences of real numbers. 

Lemma 2.1 If (aij)i^jj^j is a double indexed bounded subset of real num- 
bers, which is directed in each index, then 

SUp SUp Oij = SUp SUp Oij . 

Proposition 2.2 The partial order on partial states is directed complete. 
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Proof Given a directed subset of partial states, consider the map 

p : S(H) -> [0,1] given by p = sup i€l pi, i.e., = sup i£l pi(A) for any 

Borel subset A C R. Clearly p({0}) = and < 1. Suppose (A/)j>o 

is a sequence of mutually orthogonal subspaces. For each i € I we have by 

assumption: P;(V„>o^«) = Y, n >oPi( A n)- Le t B ij = ELoft(^n) be tne 
partial sum of the infinite series of nonnegative terms Pi(Vn>o^)- Then 
the bounded double indexed subset (■Bij)ieJ ) j>o of non-negative numbers is 
directed in each index. It follows from the lemma that we can change the 
order of taking suprema in By and thus: 

P(V„>o A «) = su Pi6/W(V„>0 A n) 
= sup ie/ supj> Bij 
= sup,-> sup ie/ By 
= SUp^Q SUp ieJ YH=oPi( A n) 
= su Pj>oJ2li=oP( A n) 
= J2 n >0P( A n) 

It follows that p is a partial state. It is easy to check that p is indeed the 
supremum of the subset (pi)igj. □ 

We denote the complete partial order of partial states (subprobability mea- 
sures) by M(S(TC)) and the set of states (probability measures) by M 1 (S(TC)). 

3 Extension of Gleason's theorem 

Consider the set D 1 of density operators, in other words, positive linear 
operators / : H — * H of trace class with trace equal to one, i.e., (x\fx) > 
for all x £ H, which we write as / > 0, and tr(/) = X^>o( e *l/ e «) = 1 f° r an 
orthonormal basis (ej)j>o of TC. In fact, the trace of a positive operator is 
independent of the orthonormal basis. Consider the linear map G : D 1 — ► 
M 1 (S(H)) given by G(f) : S(H) -► [0, 1] with G(f)(K) = tr{P K f), where 
is the orthogonal projection onto the closed subspace K G S(7i). The 
following celebrated theorem was assumed by von Neumann, then conjec- 
tured by Mackey and finally proved by Gleason. A more elementary proof 
was later provided in 0. 

Theorem 3.1 0/ The map G is onto if the dimension ofTt is greater than 
two. □ 

In order to obtain an extension of Gleason's theorem for partial states, we 
consider the notion of partial density operators. A positive linear operator 
/ : H — > H of trace class is a partial density operator if tr(/) < 1. Partial 
density operators are equipped with the usual partial ordering of operators, 
namely g < f if/ — g is a positive operator, written as: / — g > 0. Let D be 
the partial order of partial density operators. We use the following result of 
Vigier. 
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Theorem 3.2 [?[ page 51] Every norm bounded increasing filter of opera- 
tors has a supremum. 



Corollary 3.3 The partially ordered set of partial density operators is di- 
rected complete. 

Proof A directed set in D is an increasing filter with respect to the partial 
order. Since for any operator / £ D, we have ||/|| < tr(/) < 1, it follows 
that any directed set in D is bounded with respect to the operator norm 
and the result follows from the above theorem. □ 

Recall that for directed complete partially ordered sets R and S, a map- 
ping h : R — > S is an order isomorphism if h is one-to-one, onto, monotone 
(h(x) < h(y) if x < y) and reflects the partial order (x < y if h(x) < h(y)). 
We note that any order reflecting map is necessarily one-to-one: h(a) = h(b) 
implies h(a) < h{h) and h(b) < h(a), from which it follows that a < b and 
b < a, i.e., a = b. 

Proposition 3.4 // h : R — > S is an order isomorphism, then it preserves 
the supremum of directed subsets, i.e., sup ig/ /i(aj) = /i(sup ie/ aj) for any 
directed subset (aj)j g /. 

Proof By monotonicity of h we have: supj g //i(oj) < /i(supj g / Oj). Since /i 
is onto, there exists a € R such that h(a) = sup ig/ h(ai). Thus, for each 
i 6 I, we have: fa(dj) < /i(o)) < /i(supj g j aj). By order reflection, it follows 
that for i £ I, we have: en < a < supjg/Oj. Hence, a is an upper bound 
of the directed subset (aj)j g / and therefore sup ig/ a, < a, which implies 
a = supj g/ a». □ 

We will now state and prove the extension of Gleason's map to partial 
density operators and partial states. 

Theorem 3.5 For Hilbert spaces of dimension greater than two, the directed 
complete partially ordered sets of the partial density operators and the partial 
states are order isomorphic. 

Proof Consider the linear extension G : D — > AI(S(TC)) given by G{f) : 
S{H) -» [0, 1] with G(f)(K) = tr(P K f), where P K is as before the orthog- 
onal projection onto the closed subspace K € S{TL). To show that G is 
monotone, suppose g < f and assume K £ S(H). Let (ej)j>o be a complete 
orthonormal set of eigenvectors of P K . Then, 



trP^ 



(f-9) 



Ei>o^i\P K (f- 

Ei>oiP K ei\(f- 

Ep«e i =e j ( e i I (/ 



-g)ei}>0. 
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To show that G is order reflecting, assume G{g) < G(f). If, for some 
x G H with ||x|| = 1 we have {x \ (f — g)x) < 0, then by taking K to be 
the subspace generated by x and using an orthonormal basis (ej)j>o with 
eo = x, we obtain: 

G(f)-G(g) = trP K (f-g) 

= Ei>o{ei\P K (f-9)eJ 

= Ei>o(P K ei\(f-9)ei) 
= {x | (/ - g)x) < 0, 

which gives a contradiction. This proves order reflectivity, from which it 
follows that G is one-to-one. The surjectivity of G from Gleason's theorem 
by noting that for < r < 1 and all / G D, we have: G{rf) = rG(f). □ 

We note that the Scott topology of a directed complete partial order A is 
given as follows. The open subsets of the Scott topology are those subsets 
O C A such that (i) O is upward closed, i.e., x G O and x < y implies y £ O, 
and (ii) O is inaccessible by directed subsets, i.e., whenever swp ieI Xi G O 
for a directed subset (xj)j g /, then there exists i G / with xi G O. The Scott 
topology on a directed complete partial order is in general To and is the 
canonical topology with respect to which a directed set (net) converges to its 
supremum. In fact, in analogy with continuous maps between metric spaces, 
we have the following property: A mapping h : A — > B between directed 
complete partial orders is continuous with respect to the Scott topologies 
of A and B if and only if / is monotone and preserves the supremum of 
directed subsets pQ. We have thus obtained canonical topologies on D and 
M{S{TL)) with respect to which the one-to-one correspondence between the 
partial density operators and the partial states is in fact a homeomorphism: 

Corollary 3.6 The map G : D — > M{S{TL) is continuous with respect to 
the Scott topology. □ 

4 Expected values of observables 

Recall ^0] that an observable is a map r : B(R) — > S(TC)), where B(W) is 
the collection of Borel subsets of the real line, such that: 

(i) r (0) = and r(R) = H. 

(ii) r(U) and r(V) are orthogonal subspaces if U fl V = 0. 

(iii) r (Un>o U n ) = Vn>o r (U n ) for any sequence (f7 n )n>o of Borel subsets. 

The interpretation of the map r as an observable is that for the quantum 
event r{U) G S(Ti), with U C R a Borel subset, the observable takes its 
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values in U. The observable r is bounded if there exists a > such that 
r([— a, a]) = TL. The spectrum Spec(r) of r is the intersection of all closed 
subsets C C R such that r(C) = TL. By the spectral theory of operators, 
there is a one-to-one correspondence between (bounded) observables and 
(bounded) self-adjoint operators on TL pages 235, 263], which are the 
standard representation of observables in quantum mechanics. 

If p : S(TL) — > [0, 1] is a state, then the composition Q p = p a r : B(M) — > 
[0, 1] with Qp(U) = p(r(U)) is a probability distribution on R and the ex- 
pected value of r in the state p is defined by 

/oo 
tdQ r p (t), 
-oo 

when the integral exists. If r is a bounded observable and p = G(f) for 
some / G D 1 , then the expected value of r with respect to p exists and is 
given by 

E(r | p) = tr(Ar/), 

where ^4 r is the unique bounded self adjoint operator corresponding to r [101 
page 61]. 

The following question now naturally arises. Given an observable, can we 
define the expected value of a partial state, extending the classical definition 
of the expected value with respect to a state? We will show here that we can 
do this in a satisfactory way for the case of bounded observables by defining 
the expected value with respect to a partial state as an interval rather than 
a real number. 

Consider the observable r with respect to a partial state p : S(TL) — > 
[0,1]. Let m = infSpec(r) and M = supSpec(r). The composition Q r p = 
p o r : B(R) — > [0, 1] with Q r p (U) = p(r(U)) is a subprobability distribution 
on R. We put 

E (r\p)= / tdQ r p {t), 

J m 

when the integral exists. When r is bounded, we define the expected value 
of r with respect to the partial state p as the closed bounded interval: 

E(r\p)=E (r\p) + (l-Q r p (R))[m,M], (1) 

where we have used the standard notation in interval arithmetic: k+ [a, b] = 
[k + a, k + b] for a real number k and a real interval [a, b]. The second 
term in Equation ^ says that the missing probability, i.e. (1 — Q p (M)), is 
known to be distributed on [m, M] but, at this stage of computation, its 
precise distribution on this interval remains indeterminate. This uncertainty 
explains why the expected value with respect to a partial state is an interval 
rather than a real number. For a non-empty compact interval a, we use the 
notation a = [a, a] for the left and right end points of the interval. 
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Proposition 4.1 Suppose r is a bounded observable. Ifp < q, then E(r\p) 3 
E(r\q). 

Proof Let \l = Q r q — Q p . Then by assumption \i is a nonnegative measure 
with support in [m,M]. Thus: 

M 

J tdfjt(t) > /x([m,M])m, 

m 

and, hence, 

pM pM 

/ tdQ r q (t)- tdQ r p (t)>(Q g (R)-Q q (R))m. 

It follows that E(r\p) < E(r\q). Similarly, E(r\q) < E(r\p). □ 

We now aim to show that our notion of expected value has the required 
limiting properties. 

Lemma 4.2 Suppose r is any observable. If (pi)i<=i is a directed set of 
partial states with p = sup ig/ pi and k : R — > R is any measurable function 
such that k dQ r p exists, then 

/oo poo 
kdQ r p = ]im kdQ r pv 
-oo J —oo 

where the limit is understood as the limit of the net k dQ' Pi )i e i . 

Proof Let fii = Q p . and fi = Q p . Let h : R — > R be any simple function, 
i.e., h = Y^=i Xb 3 for some measurable subsets Bi C R (1 < j < n), where 
X# is the characteristic function of .B. Then, 

suPieiJ^hdm = su.p ieI J!° 00 J2 1 j=iXB j diH 

poo v-^n , 

= Lo L, i =iXB j dn 

Assume first that k is nonnegative. Since the set of simple functions h with 
< h < k is directed, from the above we obtain: 

J^ 3 ^ kdfi = sup{ J*^^ /id// | nonnegative simple h < k} 

= sup{sup ie/ h d^i | nonnegative simple h < k} 

= sup ie 7 sup{ h dfii | nonnegative simple h < k} 

= sup ie/ fc d/x; 

= lim i6/ f^kdfr. 
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This establishes the result for a positive function k. More generally, we 
split k into its positive and negative parts and put k + = max(/c, 0) and 
k~ = max(— k, 0), where denotes the constant function with value zero. 
Since J_^ o k dQ r p exists, it follows that both J_f^ o k + dQ p and A: - dQ^ 
exist as well. From k = k + — k~ , we obtain, using the result for a positive 
function: 

JZc kd V = I-oc k+ % - f-oo k ~ dfi 

= lim ie/ k + dm - lim ie / fc - 
= lim ie j k + - k~ dm 
= lim i6 7 f^ kdfXi, □ 

Corollary 4.3 If (pi)ig/ is a directed set of partial states, then for any 
observable r such that Eo(r \ sup^jPi) exists, we have: 

sup E (r | pi) = E (r \ sup^)- 
iei iel 

We can now show that the expected value of an observable with respect 
to a partial state has the desired limiting properties. Note that the collection 
of non-empty compact intervals of the real line partially ordered by reverse 
inclusion, denoted by IK, is a directed complete partial order, in which 
the supremum of a directed set is the intersection of the compact intervals 
represented by the directed set. 

Theorem 4.4 For a given bounded observable r, the expected value map 
E(r\.) : M(S{7i)) — > JM is Scott continuous. 

Proof We have already checked in Proposition 14.11 that the map E(r\.) is 
monotone. Let (pi)i£i be a directed set of partial states with p = sup iGl p^. 
Then, by Corollary 14.31 we have: 

lim ieI E{r\ Pi ) = lim ieI {E {r\ Pi )+m{l-Q r pi {R))) 

= ]mu eI E {r\ Pi ) +m(l - \im ieI Q r pi (R))) 
= E (r\p) + m(l - Q£(R))) = E(r\p) . 

Similarly, limj g / E(r\pi) = E{r\p) and thus: f] i&1 E(r\pi) = E{r\p). □ 

The representation of observables as functions of type B(R) — > S(7i) has 
serious drawbacks when studying the calculus of observables. For example, 
representing a linear combination of commuting observables and computing 
its expected value in this setting are far from straightforward; see |10( pages 
125, 163]. In fact, for this purpose, it is far simpler to represent observables 
as self adjoint operators on the Hilbert space. 
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Proposition 4.5 If A r is the self adjoint operator corresponding to the 
bounded observable r and if p = G(f) is a partial state with f a partial 
density operator, then we have: 

E (r\p) = tr(A r f). 

Proof In case, p is a state, and hence / a density operator, the result follows 
from jlUl theorem 7.24]. If p = is the trivial partial state, then / is the 
zero operator and we have Eo(r\p) = ti{A r f) = 0. Otherwise k=p(7i) > 
and p/k is a state. By linearity of G, we have: G(f /k) = p/k. Now the 
linearity of integration and the trace function implies: 

E (r\p) = kE (r\p/k) = ktv(A r f/k) = tr(A r /). □ 

Corollary 4.6 With the above notations: 

E{r\p) = tr{A r f) + (1 - tr(f))[m,M\. □ 

Since every partial density operator / corresponds to a unique partial state 
G(f), we can now write the expected value of a bounded self adjoint operator 
A with respect to a partial density operator / as 

f(A|/)=tr(A/) + (l-tr(/))[m,Af], 

where m = infSpec(^4) and M = supSpec(^4). Finally, we deduce that, 
as an interval valued map, the expected value with respect to a partial 
state depends linearly on commuting observables. For a bounded self ad- 
joint operator C, we write: mc = infSpec(C) and Mc = supSpec(C). 
Then [m^c, Mkc] = k[mc, Mc], for any real number k, and, for com- 
muting bounded observables A and B, we have: vtla+b = n^A + m B and 
M A+B = M A + M B , i.e., 

[m A+ B,M A+B ] = [m A ,M A ] + [m B ,M B ]. 



Proposition 4.7 If A and B are commuting bounded self adjoint operators, 
k and I real numbers, and f a partial density operator, then 

£(kA + lB\f) = k£(A\f) + ££(B\f). 

Proof We have: [m kA+ i B , M kA+eB ] = k[m A ,M A ] + £[m B ,M B ]. Hence, 

£(kA + lB\f) = tv((kA + £B)f) + (l-tT(f))[m kA+eB ,M kA+eB ] 

= ktr(Af) + £tv(Bf) + k(l - tt(f))[m A , M A ] + £(1 - tr(/))[m B , M B ] 
= k£(A\f)+££(B\f). □ 
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One can also compute the expected value of various functions of an ob- 
servable. For example, for a bounded self-adjoint operator A and a partial 
density operator / we have: 

S{A 2 \f) = ti(A 2 f) + (1 - tr(/))[m A2 ,M A2 ], 

where m A 2 = k 2 and M A 2 = K 2 , with k = inf{|a| : a € Spec(yl)} and 
K = sup{|a| : a <G Spec(A)}. 
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